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Abstract 

A simple relativistic quantum field model with the Yukawa-type interaction 
is considered to demonstrate that the analytic confinement of the constituent 
("quarks") and carrier ("gluons") particles explains qualitatively the basic 
dynamical properties of the spectrum of mesons considered as two-particle 
stable bound states of quarks and gluons: the quarks and gluons are con- 
fined, the glueballs represent bound states of massless gluons, the masses of 
mesons are larger than the sum of the constituent quark masses and the Regge 
trajectories of mesonic orbital excitations are almost linear. 
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I. INTRODUCTION 



The meson spectroscopy as the theory of bound states of quarks, and the phenomenology 
of the Regge trajectories (RTs) are important and interdependent subjects of investigation 
in particle physics [0-0. 

At the present time, QCD is commonly considered the true theory of the strong in- 
teraction describing all processes in the hadron world, including the mesonic spectroscopy 
0. However, being a nonlinear theory with local colour gauge symmetry, QCD is quite 
complicated from the computational point of view and corresponding methods of calcula- 
tions require great efforts in additional assumptions and ideas. In contrast to QED, simple 
and reliable methods of calculations are still missing in QCD. From our point of view, the 
main puzzle is that any acceptable and well established description of quarks and gluons 
explaining the hadronization on large distances, where the confinement of quarks and gluons 
takes place, has not been found yet. On the other hand, one may expect that a theoretical 
description of colorless hadrons considered as bound states of quarks and gluons, when the 
confinement is taken into account and the averaging over all non-observable color degrees 
of freedom is performed, can lead to a physical picture, where the quarks and gluons are 
realized in the form of some phenomenological "bricks". We suppose that a successful guess 
of the structure of these "bricks" in the confinement region can result, particularly, in a 
quahtatively correct description of the basic features of the meson spectrum. Our guess is 
that the analytic confinement realizes these "bricks". 

Following this idea, in the present paper we consider the meson spectroscopy within the 
theory of bound states of quarks and gluons. Our aim is to show within a simple relativistic 
quantum field model that the basic features of experimentally observed meson spectrum (see, 
e.g. [§]) can be qualitatively explained by the analytic confinement of quarks and gluons. 
The basic dynamic characteristics of spectrum of all mesons considered bound states of 
quarks and gluons (in contrast to the relations of the SU3 flavor symmetry) can be listed as 
follows: 

- quarks and gluons are confined; 

- glueballs are bound states of massless gluons, i.e. they are completely relativis- 
tic systems; 

- the masses of mesons are larger than the sum of masses of the constituent 
quarks; 

- the RTs of different families of mesonic orbital excitations are quite close to 
linear and their slopes practically coincide. It means that the slope is an universal 
parameter which is defined by the general nature of quark-gluon interaction. 

Obviously, these characteristics cannot be obtained in the framework of any local quan- 
tum field theory, where the constituent particles, quarks and gluons, are described by stan- 
dard Dirac and Klein-Gordon equations. From the common point of view the confinement 
plays the first role in understanding and explaining of this picture. The point is how to 
realize mathematically the conception of confinement in specific theoretical formalism? 

The standard theoretical calculations leading to linear RTs of mesonic bound states of 
quarks and gluons are based on the following assumptions. First, the quarks are accepted 
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as ordinary fermions interacting by means of gluons. Some argumentations including lattice 
calculations of the nonlinear QCD gluon dynamics are given to get a particular infrared 
behavior of the gluon propagator (D(p^) ~ for 0) which results in a linear 

increasing potential between quarks in three-dimensional space x G (see, for example, 
0]). Exactly this infrared behavior is interpreted as the quark confinement. 

Second, one of three-dimensional reductions of the relativistic Bethe-Salpeter equation is 
used. For this purpose it is necessary to overcome some mathematical problems caused by 
the singularity of the gluon kernel and an ambigously defined choice of particular reduction 
of the relativistic two-body Bethe-Salpeter equation (see, for example, As a result 

a relativized Scrodinger equation is obtained in R^-space to estimate the meson masses 
including higher orbital excitations. Hereby, the kinetic energy of the quarks looks like 
Ekin = Vp^ + cind the potential energy increases linearly. The most essential part of the 
two-body Hamiltonian takes the form 

H = \Jpi + ml + ^Jpl + ml + a\ri - rs] . 

In the present paper we demonstrate another possible mechanism explaining the above 
mentioned characteristics of mesonic spectrum, including the origin of linear RTs. In doing 
so, we use a simple relativistic quantum-field model of two scalar particles (the prototypes 
of constituent "quarks" and intermediate "gluons") with the analytic confinement. Our 
approach is based on the following assumptions: 

- the analytic confinement takes place, it means that propagators of "quarks" 
and "gluons" are entire analytic functions of the first order on the complex 
momentum p^-plane; 

- the interaction is described by the Yukawa type Lagrangian; 

- the coupling constant binding the "quarks" with "gluons" is small; 

- final bound states of "quarks" are described by the relativistic Bethe-Salpeter 
equation in one-" gluon" exchange without using any three-dimensional reduc- 
tion. 

In the present paper we shall not discuss the origin and details of the analytic confine- 
ment. Particularly, the self-dual homogeneous vacuum gluon field in QCD results in the 
analytic confinement of quarks and gluons PJ^. The model of induced quark currents based 
on this vacuum gluon field describes correctly the experimental data [p!0| . 

This paper is aimed to demonstrate the mathematical sketch of calculations of the two- 
body bound-state spectrum within the Bethe-Salpeter equation within simple QFT models, 
when the analytic confinement takes place in the weak-coupling regime. These simple rela- 
tivistic models are based on physically transparent hypotheses and can be treated by simple 
analytic methods. We claim that the analytic confinement is the basic underlying "brick" 
leading to the qualitatively correct description of main characteristics of the meson spectra. 
In any case we believe that our models represent certain theoretical interest because they 
clarify the underlying physical principles of the meson spectrum. 

For simplicity, we deal with two possible versions of the analytic confinement. In the 
first model the propagators of "quarks" and "gluons" are pure Gaussian exponents, i.e. 
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D{p^) ~ e~^^/^^, where the parameter 1/A represents the "radius" of confinement. From 
physical point of view this model is important because the eigenfunctions and eigenvalues of 
the relativistic Bethe-Salpeter equation within one-particle exchange approximation can be 
found explicitly and the obtained RTs are purely linear. One can say that this model can 
be considered a "relativistic oscillator", i.e. it is a relativistic generalization of the quantum 
nonrelativistic oscillator with equidistant spectra. 

The second model comes from a natural physical picture of standard particles with a 
relevant scale of mass m. It is supposed that the analytic confinement characterized by the 
scale A makes the propagators of particles ("quarks" and "gluons") entire analytic functions. 
As a result there are two physical parameters: the scale of confinement A and the mass of 
quark m. An interesting problem is to investigate the dependence of the meson spectrum 
on the relative parameter fi = m/A. We shall show that this model describes qualitatively 
well all above listed dynamical characteristics of meson spectra. 



II. TWO-PARTICLE BOUND STATES 

We consider the Yukawa model of two interacting scalar fields and (p{x) described 
by the following Lagrangian in the Euclidean space-time as follows 

= -$+(x)S'^(n)$(x) - ^if{x)D-\n)^{x) - ^$+(x)$(x)v9(x). (1) 

The coupling constant g is supposed to be sufficiently small. 

We postulate that the analytic confinement takes place here. It means, that propagators 
S'(p^) and -D(p^) of confined particles $ and ip are entire analytic functions in the complex 
p^-plane so that the functions S~^{z) and D^^{z) have no zeros at any finite complex z. As 
a result the equations for free fields 

S-\a)(!>{x) = 0, D'\n)^{x) = (2) 

result only in the unique trivial solutions <l>(x) = and ip{x) = 0. We call this property 
analytic confinement, i.e. the corresponding particles exist only in virtual states, so $ and 
(f are the virion fields (see 0). One can say that these fields describe constituent particles, 
e.g. $(x) and ip{x) represent scalar "quarks" and scalar "gluons", correspondingly. 

As mentioned above, we deal with two simple versions of the analytic confinement. 

First, we consider the simplest, but important case with pure Gaussian propagators 

S(p^) = ^e-ii, S(.) = ■ e-i'^''\ (3) 
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where the parameter e -C 1 implies that $-particle is much "heavier" than (/^-particle. In 
some sense this model can be considered a relativistic "oscillator" because the corresponding 
Bethe-Salpeter equation can be solved exactly and the exact solution gives linear RTs. We 
call this case the virion model. 
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The second model implies that there exists a certain physical mechanism which generates 
analytic confinement of standard particles with the initial masses m and 0. Then, their 
propagators are given in more realistic forms 



= (4) 

The confinement region is characterized by the parameter A. In the deconfinement limit 
A — > one obtains the standard propagators of massive and massless scalar particles. Within 
this model we want to analyze the infiuence of the parameter /j, on the behavior of the meson 
spectrum. We call this case the scalar confinement model. 

"Two-quark" bound states can be found in the following way. Let us consider the par- 
tition function 

III e-(*^^"*)-5(^^-'^)-^(*+*^) . (5) 

By integrating over ip we arrive at 

The term L2[$] = can be transformed as 

^2[*] = ^ JJ dXidX2 - X2)^^{X2)^{X2) 

dx jj dyidy2 ^jD{yi)J{x,yi)S{yi - y2)y/D{y2)J^{x,y2) , 
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where xi = x + ^, X2 = x — ^ and 



J{x,y) = [x + l-y^ ^ (x -^-y] = $+(a;)e2 ^$(a;), 



2V \ 2' 
J^{x,y) = J{x,-y). 

Let us introduce an entire orthonormal system {f/gd/)}, where Q = {n,l, {n}} may be 
considered as a set of radial n, orbital I and magnetic {/i} — {/ii, ...,iJ,i) quantum numbers. 
We have the conditions 

J dy UQ{y)UQ,{y) = 5qq,, Y.UQ{y)UQ'{y') = 5{y - y'). 



Q 



Then, 
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^2[*] = y E / dx Jq{x) ■ Jq{x), Jq{x) = ^+{x)VQ(dmx), 
Q 

J+{x) = Jq{x) , Vgid) =i' jdy ^D(^t/Q(y)ei^. 



By using the Gaussian functional representation we write 

2 

^ ^/dx Jq{x)-Jq{x) 
Q 

-lJ2{BQBQ)+gJ2J dx Bq{x)Jq{x) 

Q 

We substitute this representation into the partition function (1) and integrate over 
The result reads 

-\Y.<^BQBQ)-i;),\n{l-gBQVQS) 

Q 

-h E(Bq1^qq'-^^qq']Bq')+WiIB] 
U^Bqc Q'^' . (7) 

Q 

Here 
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Wi[B] = -Tr 



ln(l - gBQVQS) + ^-BqVqSBq,Vq.S 



is a functional, which describes the interactions of the fields Bq. 

The polarization operator AIIqq/ in the one-loop approximation is 

mQQ'{z)=g^Ti{VQSVQ, S) (8) 
= jj dyidy2 f/Q(?/i)An(z; yi, ?/2)f^Q'(z/2) , 

yi -y2\ a f y^- y-i 



\Ii{z- y,, y,) = ^ D{y,) S [z + S [^z - ^ D{y,) , 

where z = xi — X2 and A = {g/AnKY. Its Fourier transform reads 

Anp(!/i,!/2) (9) 
= 0^4^) ji. e'"S (. + ^) S (. - ^) 

The orthonormal system of functions {Uqly)} should diagonalize the kernel (|^). For this 
aim we consider the eigenvalue problem 

J dy' \%{y, y')UQ{y) = Egi-p^Pgiy), (10) 

where Eqi^—p^) = Eni{—p^), i.e. the eigenvalues are degenerated over the magnetic quantum 
numbers {/u}. 

If the functions {Uqiy)} are found, the polarization operator I[p{x,y) is diagonal on 
these eigenfunctions 

XIi,{y,y') = Y.EQ{-p')UQ{y)UQiy') . (11) 
Q 

We note that the diagonalization of Ilp{y, y') is nothing else but the solution of the Bethe- 
Salpeter equation in the one-boson exchange approximation. To get the standard form of 
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the Bethe-Salpeter equation, we have to introduce new functions UQ{y) = y D{y)'^Q{y) (see, 
e.g., and go to the momentum space. 
So, the partition function (|^) becomes 

U^Bqc <^ , (12) 



GQ' = GQ\-D) = l-EQ{n), p' = -D. 



The functional integral ( [T^ ) is given over a Gaussian measure defined by the operator 

We would like to stress that this representation is completely equivalent to the initial 
one d^). From physical point of view, we pass on from the world containing fields $ and 
to the world of bound states {Bq}. 

The fields {-Bq} can be interpreted as the fields of particles with quantum numbers Q 
and masses Mq, if the Green functions Gq{p'^) have simple poles in the Minkowski space 
{p^ = —Mq), in other words the equations 

1 = Eq{M'q) (13) 

define the masses Mq = Mni of two-particle bound states with quantum numbers Q = {nl) . 
These states are degenerated over magnetic quantum numbers {//}. The operator Gq^(— □) 
defines the kinetic term of the field Bq. To represent this operator in the standard form we 
expand it in the vicinity of = —Mq as follows 

1 - EQi-p') = Zq(/ + m'q) + 0[ip' + Mj)2] , Zq = -E'q{-MI) > , 

where the positive constant Zq defines the renormalization of the wave function of the field 
Bq. By using the standard renormalization relation 

Bq{p) = Zq^"Bq{p) 
we write the kinetic term in the standard form: 



(5+(p) [i-Eq{-p')\ Bq{p 
= [ip' + M'q) + 0((/ + M^q)')] BQip) 



The interaction functional in (O) reads 



Wj[gB] = Wj[g,sB] , g^^ = gZQ^'^ = ^ > 0. (14) 



It is essential, that according to ([T0| ) the effective coupling constant g^^ does not depend on 
g in any explicit way. 
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III. THE VIRTON MODEL 



Let us consider a pure Gaussian case of analytic confinement where the parame- 
ter 1/A represents the "radius" of confinement. Due to the Gaussian character of these 
propagators, the polarization kernel (|]) is quite simple 

Xfl,{y, y') = ■ e'^ ■ K{y, y') , K{y, y') = e-^(J'^-^f'+^") (15) 

with A = [g/AixK)^. The kernel K{x,y) can be explicitly diagonalized 

Kiy,y') = J2UQiy)KQUQiy'), (16) 
Q 

2n+l / 



(271 

I^Q — l^nl — 1^0 ■ \ 7= I , /^O 



A2(2 + V3), 



The eigenfunctions Ugly) are given in Appendix A. 

According to ([T3| ) the mass spectrum of two-particle bound states can be found explicitly 

Mj = = Ml + (2n + /) • 2A2 ln(2 + v^) , (17) 



Ac , /4(2 + v^)2^ 



2 



Ml = 2K^\n'^, A„ 

A V ^ / 

Thus, the Gaussian form of analytic confinement (^ leads to the purely linear and 
parallel RTs. The slope of RTs is defined only by the scale of the confinement region A 
and does not depend on the coupling constant A and other dynamic constants entering the 
interaction Lagrangian. 

Bound states exist in the weak-couphng regime A < Ac- If A <^ Ac the size of the 
confinement region is remarkably larger than the Compton length of all bound states 

^conf ~ -T > f 



A Mq 

In other words, the physical particles, described by fields Bq{x), and all physical transfor- 
mations involving them take place within the confinement region. 



IV. THE SCALAR CONFINEMENT MODEL 

The propagators of this model are given by (^. In this case the eigenvalue problem (p!0D 
cannot be solved exactly, so a variational approach will be used to evaluate approximately 
the two-particle spectra. Further we consider only the orbital excitations, i.e. the eigenvalues 
for zero radial quantum number n = 0, so that Q = {0,1, {/i}}. If {"^Q^y)} = {^«{^}(y)} is 
a set of normalized trial wave functions, then we have to calculate 

max 5^ // dyidy2 ^Q{yi)XUp{y,,y2)^Q{y2) = ei{Mf) < Ei{Mf) 

(18) 
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with p2 = -Mf. The mass of the bound state is determined by the equation (|13]): 

ei{Mf) = 1. (19) 

Note, the solution of this equation gives the upper bound to the mass Mf, because ez(M^) < 
Ei{M'^) for positive > 0. 

In the present paper we choose the normahzed trial wave function in the following form 



C, = A'+i 



(20) 



'1 + 2a)'+i 



\ 2'(/ + l)! 



M 



where a is a variational parameter. The four- dimensional spherical orthogonal harmonics 
Ti{^}{x) are defined in Appendix A. We guess that the test function in ( ^OD should be a 
sufficiently good approximation to the exact one because the kernel (0) is proportional to 



D{y) and S{y) is of the Gaussian type. 
Let us define 



^H,}ik,a) = jdx e-'^^.jD{x)^i{,}{x,a) = Ti{,}{k)I{U 



(27r)' 



where 



Here the rotational symmetry D = D{y'^) = D(Y'^) has been used. One obtains 



24+31 

4 



(21) 



Mo 



Substituting representations 
arrive at 

dk 



A2(l + a) ■ 

]) and (l2l|) into (|18D and after some calculations we 



<^ J (27r) 



{m} 



ei{Mf) = / max /-^ ^ <l,|^}(fc, a) G (k + ^) G (k - ^) ^i{,}{k, a 



P 



2A 



■ max <j [4c(l - c)f+' jjdtds e-(*+^)(A' R^it, s,x) \ = 1- 



(22) 



Here 



Riit, s, x) = JJdudv e-" ^ {uwYFiib, x), 





-- 




I db 
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FIGURES 

85.00 q 
20.00 z 
15.00 - 
10.00 - 
5.00 - 

0.00 n 

0.0 3.0 -1.0 6.0 B.O 10.0 12.0 

Fig. 1 

The spectrum of two-particle bound state (Mq/A) at zero orbital 
angular quantum number / = and A = 0.01 as function of the 
mass n = m/A of the "constituent" particle. The dashed curve 
depicts the asymptotical line 2/1. 




with p 



[iM, 0, 0, 0) and 

/ n \2 



A 
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-M2, b 

m 



t + s + 2c{u + w) 
2A' 



K 



and 



V47rA. 

The mass M; of the two-particle bound state in the one-parameter variational approxima- 



tion should be determined from the variational equation (p2[) . We have solved this equation 
numerically for different values of parameters A, / and ^. The obtained results are plotted 
in Figs. 1-3. 




Fig. 2 

The spectrum of two-particle bound state as function of the or- 
bital angular quantum number / for several values of /i = m/A = 
{0.05; 0.2; 1.0; 5.0} at fixed effective coupling constant A = 0.01. 
For convenience the difference {Mi/Af - {Mo/Af is plotted. 
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Fig. 3 

The spectrum of two-particle bound state as function of the 
orbital angular quantum number / for several values of A = 
{0.001; 0.01; 0.1} at fixed fi = 1.0. For convenience the differ- 
ence {Mi/Af - {Mo/Af is plotted. 



V. DISCUSSION 

Now let us make sure that the Scalar Confinement Model correctly describes the basic 
features of the meson spectrum, when mesons are considered bound states of quarks and 
gluons. In Fig.l we show the mass Mq/A of the lowest bound state Bq with quantum 
numbers Q = (0, 0) = as the function of the mass of constituent quarks n = m/A for small 
coupling constant A = 0.01. 

First of all, one can conclude that there exists a bound state in the case /i = 0, i.e. with 
m = 0. In other words, provided that the analytic confinement takes place, two massless 
particles can be coupled into a stable bound state (see also ||12|). This situation cannot be 



realized under any circumstances in the local quantum field theory. 

Second, the mass Mq > 2m for ^ < 50 ~ 80, i.e. it exceeds the sum of masses of 
two constituent quarks. Besides, these bound states are stable. It means that these heavy 
mesons are relativistic systems and therefore, the nonrelativistic Schrodinger equation is not 
adequate tool to describe these heavy mesons. Probably the masses of real heavy quarks (see 
Table 1) are not sufficiently "heavy" to use nonrelativistic description. The conventional 
regime Mq < 2m appears for very large yU > 50 ^ 80. In the deconfinement limit A ^ 
our rough variational estimation ( p2D gives the qualitatively correct behavior (for details see 
Appendix B) 



Mo = 2m - (-^y 16nmC^, 
KAnmJ 



i.e. we get the standard nonrelativistic (the coupling constant g/iirm is small!) bound state 
under the Coulomb potential. 
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Quarks 


Mass (Mev) 


Meson 


Mass (Mev) 


u 


1^5 


p = uu 


770 


d 


3^9 






s 


75H-170 


K* = ss 


900 


c 


1150^1350 


Vc = cc 


2980 


b 


4000^4400 


'y = bb 


9460 



Table 1. 

The experimentally observed masses of light and heavy quarks q 
and mesons qq . 



In Table 1 we list the masses of quarks and mesons as bound states of corresponding 
quarks. One can see that the meson masses are larger than two masses of the constituent 
quarks. The Scalar Confinement Model discribes this property of meson spectrum. 

In Figures 2 and 3 we plot the Regge trajectories, i.e the squared masses of orbital 
excitations at different values of the small coupling constant A and masses of constituent 
quarks fi. One can see that the dependence of RTs on the parameters A and fi is not 
strong. Besides, the RTs are not linear, although the linearity begins in fact for / > 4 6. 
The asymptotical behavior of the Regge spectrum for sufficiently large / can be obtained 
analytically and coincides with the exact solution of the Virton Model (p!7|): 



Mf ~ / ■ A^2 ln(2 + V3) for / ^ oo. 



(23) 



On the other hand, the dependence on / is quite smooth, so that three or four lowest 
orbital excitations can be approximatly considered almost linear. Therefore, they can be 
well approximated by linear trajectories on a short interval / = -^ 3. 

A recent analysis of the experimental data shows (see [§]) that the RTs of different 
meson and baryon families are approximately linear and their slopes deviate slightly around 
a constant value, although the quark configurations and quantum numbers of these hadronic 
families are considerably different. Note, the analysed data in ^ have been obtained for 
low orbital momenta / = ^ 3 only. 

Nevertheless, one can conclude that the slope of RTs weakly depends on specific details 
of hadron internal dynamics and may be considered an universal characteristic which is 
dictated only by the general properties of quark-gluon interactions. 

Precisely this qualitative picture takes place for our models with analytic confinement. 
Thus, we have sufficient ground to claim that the analytic confinement realizes these general 
properties and leads to the approximate linearity of RTs for meson families. 

Therefore, the Scalar Confinement Model should qualitatively describe the correct be- 
havior of orbital excitations. We can compare the experimental data with our calculations. 
In the RTs of tt-, K-, p-, uj-, (p- and i^'*-meson families have been fitted on the interval 
< / < 3 by curves = uiq + mil + with appropriate constants rrij. We have 
also performed fittings of our data obtained for / = 1,2,3. By comparing corresponding 
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coefficients of fittings for vr- and i^-meson trajectories with our fit coefficients obtained at 
different values of parameters /i and A we have found the possible values of the confinement 
scale A. The result is given in Table 2. One can see that A depends weakly on /i and A in 
wide ranges. Our semiquantative result indicates that the confinement scale parameter is 
around Kconf ~ 500 MeV . 





A 


A^ (MeV) 


Ak (MeV) 


0.2 


0.01 


562 


504 


1.0 


0.001 


558 


502 


1.0 


0.01 


539 


487 


1.0 


0.1 


502 


452 


5.0 


0.01 


506 


456 



Table 2. 

The scale of confinement A calculated from the comparison with 
real spectra of tt- and i^'-meson families given in |^. 
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Appendix A 



Consider the kernel 

K = K{x, y) = e-'^^'+2''^2'-"?/' , a>b (24) 

with 

2 

TtK = Jdy K{y, y) = j dy e-^'~^~''^y' = -^^—g^, < oo . 

The eigenvalues with quantum numbers Q = {nl{^}} = {nl{fii...fii}} and eigenfunctions of 
the problem 



dy K{x,y)UQ{y) = kqUq{x) 
can be solved explicitly. The eigenvalues are 



\ 2n+l 9 
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The eigenf unctions are 

Uq = Urui,} (y) = N^i Ti{,} {y) Lf^'^ {2(3y^) e"^^' . (26) 
Here L'l+^\x) are the Laguerre polynomials and 



, J2HI + 1) ^ , 



Tin + l] 



IT \T{n + l + 2) 

The functions 

Ti{f,}{y) = Ti{^}{ny)\y\^, Uy = ^, \y\ = ^ 



U{^,}yy) = ^i{t,}y'^y)\y\ , ny 
satisfy the conditions 

ET,{^}K)T,{,,}(n2) = ^Cliin.n,)), =/ + l, 

where C/(t) are the Gegenbauer polynomials and 

Besides, the following relation takes place 

I A T,|,}(y)F(y2)e-'=^ = (-Z)' T,{,}(fc)J(fc2), (27) 



,4+2« 



Appendix B 



Consider the variational problem (p2D in the deconfinement limit A ^ for / = 



Aomax^c(l-c) rfttis e 1^)^*+^) (28) 



^ [t + s + 2c(M + t;)]2 

Here Mq is the mass of the lowest bound state and the effective coupling constant is given 
as 

(47rm) ^ 
13 



Going to the new variables 



A2 A2 A2 



we rewrite ([28| ) 



4Ao max < 



„2 



4m^ 



^ 

dudv 



M 



Obviously, the limit A ^ exists if Mq < 2m. Besides, since Ao < 1 and 1 - ^ < 1, the 



main contribution into the integral over dx comes from large x, so that the inner integral over 
dy can be explicitly taken on the extended interval {— cxd, cxd} without any loss of accuracy. 
Thus, we get 



Mo V 1 - MoV4m2 " ' ' ^^^^ 

C = max J ^ [dx e"^' [[- ^"^^ ,,,,, I = 0.31923.... 

o<5<oo p i + 2^{u + i;)]3/2 J 

By solving (pOf ) we obtain the mass of the lowest two-particle bound state in the deconfine- 
ment limit A — > as follows 

Mo = 2m - levrmC^ + O(A^) . (30) 

This behavior is in accordance with the conventional nonrelativistic picture of the attractive 
Coulomb interaction when the binding energy of two-particle bound system is negative. 
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